SPO11-C631T Gene Polymorphism: Association With Male Infertility and an in Silico-Analysis

Objective: To investigate the association of C631T single nucleotide polymorphisms in SPO11 gene with male infertilityfollowed by an in silico approach. SPO11 is a gene involved in meiosis and spermatogenesis process, which in humans, this gene is located on chromosome 20 (20q13.2-13.3) with 13 exons. Materials and methods: In a case-control study, 200 blood samples were collected from the IVF center (Kashan, Iran) including; 100 infertile and 100 healthy control men. SPO11-C631T were genotyped using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method.The effects of C631T transition on the structure of mRNA and protein of SPO11 was evaluated by bioinformatics tools. Results: Our data revealed that all subjects were wild-type homozygous inC631T positionsand just a sample from fertile group was heterozygousin C631T (OR: 0.3300, 95% CI: 0.0133 to 8.1992, p = 0.4988).Our in silico-analysis revealed that C631T transition could make fundamental changes in the structure of the mRNA (Score: 0.1983) and protein (PROVEAN Score: -3.371; Reliability Index: 4; Expected Accuracy: 82%) of SPO11. Also, C631T substitution could change the aggregation prone regions of the SPO11 protein (dTANGO = 209.99). Conclusion: So even though the SPO11-C631T don’t increase the risk of male infertility, it could be deleterious for themRNA and protein.


Original Article
A possible candidate for genetic susceptibility to spermatogenic failure is the SPO11gene.The human SPO11 gene consists of 13 exons and is located on 20q13.2-13.3(Figure 1A) (12).This gene encodes a critical protein with 396 amino acids.Two major SPO11 isoforms are generated by alternative splicing: SPO11α (exon 2 skipped) and SPO11β.SPO11 is a critical protein in meiosis process.Recombination during meiosis increases genetic diversity, therefore the central purpose of recombination is to promote accurate segregation of homologous chromosomes to generate haploidcells (13).During meiosis, programmed double-strand breaks (DSBs) occur, which leads to interhomolog recombination.Meiotic DSB repair is remarkably robust, which exist in every cell to repair spontaneous DNA damage.Meiotic recombination mechanisms have been deducedprimarily from studies in budding yeast (14)(15)(16).DSBs are performed by the SPO11molecule, and thenthe molecule is removed from DNA, the 5'ends are cleaved to generate 3' single-stranded tails (17).SPO11acts as a dimer, and one molecule remains covalently attached to each DSB end upongeneration of the break (18).The catalytic activity for DSBs formation seems that SPO11 presumably acts as a transesterase rather than as an endonuclease (19).There is a tyrosine residue at active site of the protein that has a crucial role in DSB (18).Yeast SPO11 mutants show a range of phenotypes from partial loss of function, to complete loss of DSB formation.Varying levels of synaptonemal complex (SC) defects are also observed, highlighting the importance of DSB formation in homologous chromosome synapsis.A previous study investigated the association of five single nucleotide polymorphisms (SNPs) in SPO11 gene with male infertility (20).These five SNPs included rs28368062, rs28368064, rs79564060, rs23736832 and rs28368082.The rs28368082 (C631T) SNP leads to Arginine to Tryptophan substitution at codon 211 (Arg211Trp) located on exon7 (Figure 1A).In the present study, we investigated the association of C631T substitution in SPO11 gene with male infertility followed by a novel in silico-analysis.

Materials and methods
In a case-control study, all subjects were randomly selected from the IVF center (Kashan, Iran).The infertile patients were selected from couples attending the infertility clinic who had a history of infertility above 1 year.After signing an informed consent, a detailed medical and reproductive history was obtained from all subjects, including reproductive history and infertility evaluation of the female partner.Samples were first screened for genetic and familial diseases.100 patients were finally selected as idiopathic infertile case group.100 healthy fertile male volunteers served as the control group.The semen analysis for sperm concentration, morphology, and motility performed following the World Health Organization (WHO) criteria (21).Finally 2 ml blood was taken from each individual and preserved at -20 °C.All the participants' informed written consent and the study was approved by the Medical Research Ethics Committee of the Kashan University of Medical Sciences on 30 August 2011 (reference number 9254).

Analysis of C631TSNP of SPO11
Human genome was extracted and showed a major band with low mobility on agarose gel.However, using the extracted genome as template, SPO11 gene containing C631T single nucleotide polymorphism with the approximate size 400 bp was amplified by utilizing the conditions which were described in materials and methods section.EcoRI cuts only the C allele, which resulted in 245 bp and 155 bp fragments.All PCR products of SPO11 fragment (400 bp) containing C631T were digested by EcoRI.The absence of polymorphism in C631T causes full enzymatic digestion of amplified fragment whereas heterozygote samples produced 3 bands on agarose gel (Figure 1B).In 199 samples in the present study, EcoRI produced two fragments with aforementioned lengths.Only one sample from control group was heterozygote in C631T (OR: 0.3300, 95% CI: 0.0133 to 8.1992, p = 0.4988).

Structural analysis
The physicochemical properties of SPO11 protein, determined using ProtParam, indicated that the SPO11 protein contains 396 amino acids.The normal known SPO11 protein has a predicted molecular formula of C 1988 H 3195 N 539 O 576 S 21 , molecular mass of 44536.7 Da and an isoelectric point (pI) of 9.05.In the normal protein, the total number of negatively charged residues (Asp and Glu) is 41 and the total number of positively charged residues (Arg and Lys) is 50.The aliphatic index of the SPO11protein had been calculated to be 96.77.The instability index of the SPO11 protein in normal was determined to be 59.99.The Grand average of hydropathicity (GRAVY) was calculated -0.088 for normal SPO11.The estimated half-life of the SPO11proteininnormal is 30h (mammalian reticulocytes, in vitro), (Expasy/ProtParam server).However these parameters for 631T phenotype are summarized in table 1.
The secondary structure of the SPO11 protein was evaluated by the Bioinf server (Figure2).Our data revealed that the secondary structure of the SPO11 protein differs between 631C and 631T phenotypes, especially in the helix, strand and coil motifs.Analysis of three-dimensional structure of SPO11 protein by Accelrys DS Visualiser 4.0 revealed that residue 211 located in a beta-sheet motif at the surface of SPO11 (Figure 3).Additionally, the data obtained by a Ramachandran plot confirmed the structural stability for both 631C and 631T phenotypes of the protein (Figure 4A).Also data revealed that the peak of the hydrophobicity for 631Cwas -0.02 whereas this index shift to 0.7 for 631T (29) (Figure4B).Analysis of C631T transition bySNP effect server indicated that C631T does not affect the chaperone binding tendency (dLIMBO: 0.00) and amyloid-forming regions (dWALTZ: -27.61) whereas the substitution increases aggregation prone regions in SPO11 protein (dTANGO: 1107.96).Data from SNAP server indicated that the C631T substitution may have damaging effects in SPO11 structure (Prediction: Nonneutral, Reliability Index: 3, Expected Accuracy: 78%) consistent the data from SIFT server (PROVEAN score: -3.371, Prediction (cutoff = -2.5):Deleterious).While the data from PolyPhen2 showed that C631Thave a low probability of damaging (Score: 0.014, Sensitivity: 0.96, Specificity: 0.79).Predicting C631T effects on local RNA secondary structure of SPO11 revealed that the SNP make fundamental changes on the secondary structure of mRNA (Distance: 0.1080, P-value: 0.1983; the P-value less than 0.2 is significant structural change) (Figure 5).

Discussion Discussion Discussion
SPO11 is a critical gene and is involved in the recombination process of genes in meiosis (30,31).Although many studies have been done on the biological role of SPO11 gene but its role in infertility has not been studied extensively yet.A common single nucleotide polymorphisminSPO11 changes a cytosine residue to thymine (631C→T) which leads to Arginine to Tryptophan substitution at codon 211 in SPO11 protein.In this study, firstly we investigate the association of 631C→T with male infertility in an Iranian population.Our results showed that only one of the fertile subjects had CT genotype in C631T location and the rest of the subjects (100 infertile and 99 fertile men) had CC genotype.In 2005 some SNPs was studied in 48 male azoospermia and 48 oligospermia in USA.Results showed that one of the azoospermisubject had a glutamate changed to lysine (Glu9Lys) in exon one another had an aspartate changed to glutamate (Asp277Glu) (32).In 2011 the relation of five different SNPs (rs28368062, rs28368064, rs79564060, rs2373683 and rs28368082) with infertility was studied in Chinese men.Only one of the SNPs (rs28368082) was related to infertility and the incidence of CT genotype in rs28368082 location in infertile men was higher than fertile men (p = 0.023).Also the frequency of T allele in rs28368082 location was higher in infertile than fertile men as well (p = 0.025) (20).In 2014 Ghalkhani et al. performed another study in Iranian population.They recruited a total of 113 infertile men, including: 58 azoospermic men and 55 oligozoospermic men.They selected the subjects from three regions; including: North (19 azoospermia and 17 oligozoospermia), Center (20 azoospermia and 20 oligozoospermia), and South (19 azoospermia and 18 oligozoospermia) of Iran.They reported that C361T was significantly associated with male infertility in all populations except oligozoospermic cases from the Center region (33).The inconsistency between the studies could arise from to race, geography or genetic differences of the study population.
In meiosis, recombination of genes occurs when homolog chromosomes separate from each other to produce haploid gametes (34,35).The expression of SPO11 gene increases during meiosis (36).In yeast, by breaking the double strand DNA, SPO11 facilitate the process of recombination.It been shown that in rats, a disrupted SPO11 gene can causes male infertility (37).Any disturbance in the recombination process of meiosis counteracts the proper dissociation of the homolog chromosomes and leads to aneuploidy and making defected gametes.Therefore any change or mutation in SPO11 gene can cause infertility in men (17,38).
Effects of missense mutations on protein features usinginsilicoprediction can be determined even without a well-defined protein structure (39).Study the secondary structure of SPO11 showed that Arg211Trp substitution makes some changes in the structure of the coils, strands, and helixes in the protein (Figure 2).The physicochemical properties of SPO11 protein, detected by ProtParam server, showed that the SPO11 protein contains 396 residues.The wildphenotype of SPO11 has a predicted molecular formula of C 1988 H 3195 N 539 O 576 S 21 , molecular mass of 44536.7 Da and an isoelectric point (pI) of 9.05; the molecular mass of the 631T phenotype increases to 44566.7 Da.Because we found that Arg had five carbon atom less than Trp, the molecular weight of the normal SPO11 protein is lower than that of the mutant protein resulting in an increase in the ratio of pI to MW for mutant protein.
In the normal protein, the total number of negatively charged residues (Asp and Glu) is 41 and the total number of positively charged residues (Arg and Lys) is 50 whereas these ratios in mutant protein are 41 and 49 respectively.Therefore pI shifts from 9.05 for normal protein to 8.97 for mutant SPO11.The estimated half-life of the SPO11 proteinin normal is 30 h (mammalian reticulocytes, in vitro).The hydrophobic pattern, determined using the Kyte and Doolittle hydrophobicity scale (29), showed that the average hydrophobicity of Arg211 in the normal protein equals -0.02 whereas this index shifts to 0.7 for Trp211 in the mutant protein (Figure 4A).Therefor the Trp211 as a hydrophobic residue tends to cluster in the interior of the molecule that this change may be deleterious for the protein (40).SIFT operates by utilizing both physical properties of amino acids and sequence homology of varied orthologus.The function of PolyPhenis based on physical and comparative approaches to evaluate the impact of an amino acid substitution on the function and structure of a protein.SNAP prediction is based on evolutionary information, biophysical properties of substitution and structure of protein and annotations (39,41).The data from SNAP and SIFT showed that C631T might be deleterious for function and structure of SPO11.Whereas PolyPhen predicted the C631T transition as a benign substitution analysis of C631T substitution by SNPeffect server indicated that the substitution increases aggregation prone regions in SPO11 protein (dTANGO: 1107.96).The total TANGO score for the protein is 1107.96.Mutations can increase (dTANGO > 50), decrease (dTANGO < -50) or not affect aggregation propensity (dTANGO between -50 and 50).For C6311 substitution, dTANGO equals 209.99 which means that the mutation increases the aggregation tendency of SPO11 protein.mRNAscontaining different nucleotide at SNP position may vary in their interactions with cellular molecules involved in mRNA transport, maturation, translation, or degradation (42).Analysis the effects of C631T transition on secondary structure of mRNA showed that the SNP could alter the structure and minimum free energy of mRNA.Therefore C631T might alter the SPO11 gene expression.
Due to the effects of gene-gene interactions and environmental factors, we need a larger group of case and control in order to obtain accurate data.Although our in silico-analysis lead to better understanding of the effects of C631T on the SPO11, but more studies are essential to evaluate the influence of the environmental factor on sperm parameters of men with different SPO11 variants.

Conclusion Conclusion Conclusion Conclusion
According to our research, although C631T polymorphism in the SPO11 gene is not associated with idiopathic male infertility in the Iran; Kashan population, but in silico-analysis of SPO11 indicated that the SNP could make fundamental changes in the protein and mRNA structure of SPO11.Therefore further studies with larger sample size and diverse ethnic populations are required to confirm the general validity of our findings.Since we only studies one SNP, further studies are required to investigate the possible effects of other mutations (alone or in combination) on SPO11 function.Also functional studies on these polymorphisms are needed to clarify the molecular mechanisms of the observed results.

Figure 1 :
Figure 1: Figure 1: Figure 1: Figure 1: Location and genotyping of C631T.A) Human SPO11 gene map was deduced from NCBI data bank, which contain 13 exon, the location of C631T SNP on exon 7, was shown by arrowhead; B) Schematic and results of the PCR-RFLP.Digestion map for 400-bp PCR fragment by EcoRI enzyme [lane M = DNA marker; lane1 = PCR product as a control; lane2 = the sample with CT genotype containing 400-bp, 245-bp and 155bp bands; lane3 = the sample with CC genotype containing 245-bp and 255-bp bands].

Figure 3 : 3 : Figure 3 :Figure 3 :
Figure 3: Figure 3: Figure 3: Figure 3: DNA binding site of SPO11 protein.Residue 211 and ligand were shown by arrowheads.The residue 211 located in a beta-sheet motif on the surface of the protein.